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L  NUMERICAL  MODELLING  OF  CLIMATIC  WATER 
MASS  DYNAMICS  IN  THE  BLACK  SEA 

R  is  a  pity  diat  the  experimental  investigadons  of  the  cunents  in  the  Black  Sea  are  limited  in 
scope  and  too  non-systematic.  In  general,  we  at  present  know  only  the  basic  rules  governing  the 
large-scale  circidadonH*  Numerical  nx^elling  is  therefore  the  only  possible  way  to  reconstruct 
the  ste^y  circulation  features  in  the  Black  Sea.  The  main  problem  solved  so  far,  using  diagnostic 
and  ad^tive  "vxirfg  of  the  Black  Sea,  is  the  calculation  of  cunents  from  experimental  observations 
of  the  hydiologicd  fields.  These  models  can  be  inteipreted  as  methods  for  hydrodynamic  analysis 
of  tte  eaqjerhnental  observatons. 

The  dynamic  method  (a  simple  geostrophic  model)  developed  by  Sandstrom  and  Halland- 
Hansenl^l  at  die  start  of  this  centmy  has  long  b^n  the  only  practical  method  for  calculation  of  die 
currents.  This  method  has  some  important  shortcomings,  but  due  to  its  simple  in^lementation  is 
widely  used  in  oceanolo^cal  worlc  A.  Sarkisyan  developed  the  first  diagnostic  model  in  the 
sixties,  used  to  perform  c&gnostic  calculations  of  the  large-scale  currents  in  the  North  Atlantic. 
Numerical  modelling  of  water  mass  dynamics  in  the  Black  Sea  was  first  performed  in  the  former 
USSR  at  the  begimiing  of  the  seventies.  Since  then  numerous  diagnostic  (and  prognostic) 
investigations  of  die  currents  in  this  basin  have  been  p^ormed.  New  and  more  accurate  models 
have  been  developol,  and  the  methods  of  the  diagnostic  calculations  have  been  inqiroved.  The 
main  results  of  th^  investigations  of  the  Black  Sea  are  published  in  various  monographsl^* 
Some  results,  obtained  at  the  Computational  Center  of  the  Siberian  Branch  of  ^e  USSR 
Academy  of  Sciraces,  under  the  guidance  of  academician  G.  Marchuk,  are  published  in^^l  as  well. 

The  first  step  in  the  dia^ostic  investigations  of  the  Black  Sea  has  included  diagnostic 
current  calculations,  using  experimental  data  on  sea  water  density  distribution,  and  based  on  simple 
modsls  and  ^arse  grids.  These  investigations  have  proven  the  realistic  character  of  the  dynamic 
qiproach  (e^edaUy  if  the  dynamic  heights  are  referenced  to  the  bottom).  They  have  tiiown  that 
tte  homogeneous  ocean  models  are  groundless  and  have  confirmed  the  main  features  of  the  large- 
scale  water  circulation  in  the  sea,  obtained  using  the  dyiuunic  method.  At  the  same  time  however, 
these  investigations  have  demonstrated  some  important  shortcomings  of  the  dynamic  method. 
A^tiomdly,  die  calculations  have  provided  the  first  substantiated  3D  fields  of  vertical  velocityl^l. 

The  next  logical  step  has  been  to  perform  diagnostic  calculations  of  the  currents,  using 
small-interval  gri&  and  higher-level  models.  The  1**  or  0.5^  grids,  used  earlier  for  the  diagnostic 
calculations  of  Ae  currents  in  the  Black  Sea,  make  it  possible  to  describe  (with  significant 
approximations)  t^y  the  most  general  features  of  the  circulation.  At  these  sparse  grids  even  the 
most  important  dynamic  features  -  including  the  General  Black  Sea  Current  (GBSC)  -  are  often 
lost  The  various  simplifying  assumptions  used  in  the  earlier  models  have  been  the  cause  of 
significant  inaccuracies.  To  overcome  these  shortcomings,  new  and  more  accurate  diagnostic 
r^els  have  been  developed  and  tested,  and  the  methods  of  performing  diagnostic  calculations 
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have  been  improved.  The  high  resolution  grids  and  more  accurate  models  used  allow  study  of  not 
only  die  dicidation  of  the  sea  as  a  whole,  but  also  investigations  of  smaller  regions,  characterized 
by  complex  dynamic  processes  and  strong  spatial  variability  of  the  currents. 

Use  of  high  resolution  grids  is  a  common  trend  for  all  numerical  models  of  ocean 
dynamics.  Besides  the  evident  advantages,  the  increased  resolution  causes  some  new  problems  to 
appear.  For  instance,  when  diagnostic  c^culations  are  performed  on  small-interval  grids,  the 
quality  of  the  input  data  is  an  inqxirmt  problem.  Significant  errors  are  caused  by  die  insufficient 
statistical  reliability  of  the  climatic  ^stributions  of  temperature  and  salinity,  the  lack  of 
synchronization  in  the  "quasi-synchronous"  hydrological  surveys,  the  experiment^  measurement 
errors,  and  the  inconsistency  of  the  input  fields  of  density,  winds  and  bottom  relief.  It's  clear  that 
die  use  of  the  traditional  data  filtration  procedures  can  not  solve  these  problems,  because: 

(a)  filtration  smooths  down  both  the  high-fiequency  noise,  and  the  real  large-scale 
features  of  the  fields,  e.g.  streams; 

(b)  the  numerous  types  of  filters  used,  although  always  claimed  to  be  unbiasing,  are  in 
fact  a  source  of  subjective  and  formal  data  modifications,  not  reflecting  the  thermo-  and 
hydrodynamical  features  of  the  real  fields; 

(c)  die  problem  of  hydrodynamical  adjustment  of  the  fields  isn't  being  solved  at  alL 

The  "defects"  of  the  input  data  and  the  shortcomings  of  the  methods  used  for  processing 
hydrological  information  have  stimulated  the  development  of  new,  more  general  hydrodynamic 
methods  of  experimental  data  filtration  and  adjustmentl^i.  For  die  Black  Sea  these  investigations 
have  been  initiated  by  Sarkisyan,  Demin  and  Trukhchevn^l  and  have  been  performed  under  the 
fiitei^ovemmental  progtam  "Sections"  and  the  topic  "World  Ocean"f^»  5. 6. 1 IJ. 

The  main  feature  of  this  new  approach  is  that  it  cmnbines  the  most  important  practical 
advantages  of  the  diagnostic  models  (being  of  practical  significance  due  to  the  use  of  experimental 
data)  and  of  the  prognostic  models  (p^ucing  hydrodynamical  adjustment  and  filtered  fields).  At 
the  same  time  the  new  approach  doesn't  have  die  shortcomings  of  the  older  models  -  inconsistency 
of  the  fields  and  the  presence  of  noise  in  them  (for  the  diagnostic  models),  and  idealization  and 
smoodiing  down  of  the  large  gradients  (for  the  prognostic  niodels).  The  new  approach  is  based  on 
die  successive  triplication  oi  both  diagnostic  and  prognostic  models.  The  temperature  and  salinity 
fields,  obtained  ^m  the  observations,  and  the  currents,  obtained  from  the  diagnostic  models,  are 
used  as  the  initial  approximation  of  the  prognostic  model  The  time  integration  process  in  die  latter 
isn't  performed  until  a  fully  stationary  solution  is  obtained,  as  is  usuaUy  done  in  the  prognostic 
models.  The  integration  is  performed  only  until  the  completion  of  the  b^c  reconstruction  of  the 
fields.  This  process  takes  place  at  the  adjustment  stage  of  the  integration  (semidiagnostic 
calculations)  due  to  die  high-speed  wave  processes,  and  results  in  hydrodynamic  noise  filtration. 
If  die  integr^cm  process  is  continued  longer,  die  solution  gets  worse,  its  correlation  with  the  initial 
data  is  lost,  and  the  final  result  is  close  to  that  of  the  ordinary  prognostic  models.  This  is  the 
reason  for  limiting  the  integration  time  in  the  adaptive  method.  The  "worsening"  of  the  solution 
can  be  expressed  by  excessive  smoothing  of  the  structures,  and  even  loss  of  some  qualitative 
features,  present  in  &e  input  data.  The  adaptation  process  is  monitored  through  the  kmetic  energy 
variations  versus  time,  and  by  test  calculations.  The  time  dependence  of  kinetic  energy  shows  the 
transition  from  the  stage  of  fast  distribution  adjustment  due  to  wave  processes  (expressed  by 
oscillations  of  kinetic  energy)  to  the  stage  of  slow,  smooth  changes  (Figuro  1). 

This  problem  is  similar  to  the  initiation  problem  in  meteorology,  but  it  has  some  additional 
complexity  due  to  the  fact  that  one  of  the  fields  (the  currents)  isn't  measured  and  has  to  be 
reconstructed  by  calculations.  The  initiation  of  the  hydrological  fields  requires  reconstruction  of 
Ae  current  distribution  from  the  given  fields  of  temperature  and  salinity,  with  some  adaptation  of 
all  three  fields  to  the  equations  and  boundary  conditions  used  by  the  model,  and  to  the  shape  of  the 
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sea  basin.  The  experience  gained  up  to  now  shows  that  the  typical  model  time  for  adaptation  for 
the  Black  Sea  is  approximately  seve^  days.  The  main  adaptation  mechanism  is  by  inertial-gravity 
waves.  As  a  result  of  this  adaptation,  "instantaneous"  or  "tozen"  hydrological  fields  are  obtained, 
corresponding  to  the  tempor^  characteristics  of  the  input  temperature  and  salinity  information 
(annud,  seasonal,  mondily,  or  simply  reflecting  one  specific  condition  in  time). 

Referencesf^*  5, 6]  contain  the  results  of  numerical  experiments  on  the  adaptation  of  main 
hydrologic^  fields  at  the  western  part  of  die  Black  Sea.  They  have  been  performed  by  3D-models 
of  varying  complexity,  using  grids  with  a  horizontal  step  of  5  miles.  Several  interesting  features  of 
the  eddy  circulation  in  this  region  have  been  obtained,  describing  its  synoptic  character  and  the 
seasonal  variations.  In  the  present  work  we  discuss  only  the  results,  obtained  by  a  model,  using 
the  full  equations  of  ocean  thermo-  and  hydrodynamics,  and  intended  for  adaptation  of 
temperature,  salinity  (i.e.,  density)  and  current  fields.  This  model  has  been  used  earlier  for 
calculations  at  the  western  regions  of  the  seal^*  ^^1.  The  aim  of  the  investigations  is  to  build  the 
climatic  (seasonal  and  annual)  distributions  of  these  fields  for  the  whole  sea  basin.  The  horizontal 
grid  steps  are  0.25®  and  0.5®  along  the  parallels  and  the  meridians  respectively.  Twenty  five  levels 
(0,  5,  10,  15,  20,  25,  30,  50,  75,  lOO,  125,  150,  200,  250,  300,  400,  500,  600,  800,  1000, 
1200,  1400,  1600,  1800,  2000  m)  have  been  investigated  in  the  calculations.  The  input 
temperature  and  salinity  fields  are  built  using  about  50,000  observations,  covering  a  period  of 
iq)proximately  100  years.  This  information  has  been  prepared  by  Dr.  G.  S.  Dvoryaninov  and  his 
colleagues  from  the  Marine  Hydrophysical  Institute  in  Sevastopol— the  T-S  data  have  been 
seasonally  and  annually  averag^  at  the  centers  of  rectangular  cells  of  0.5®  by  0.5®.  The  spring 
season  includes  March,  April  and  May,  the  summer  season  -  June,  July  and  August,  etc.  To 
improve  the  resolution  we  have  additionally  interpolated  the  source  fields  along  the  parallels.  No 
preliminaiy  ^tration  of  the  data  has  been  perfom^ 

The  numerical  model  is  based  on  the  complete  system  of  nonsimplified  geophysical 
hydrodynamics  with  the  approximations  commonly  used  in  oceanology  (in  p-plane): 
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AU  designations  sre  generally  Bcceptedi^* 

There  is  one  important  feature  of  this  model.  In  contrast  to  other  numerical  models  of 
ocean  circulation  (where  the  total  mass  transport  function  is  used  as  the  integral  function),  in  this 
case  an  equation  of  si^ace  level  is  consider^  A  similar  integral  function  .is  applied  rarely  as  a 
rule,  in  diagnostic  calculations.  The  approximation  of  the  input  equations  by  Fiadeiro  and  Veronis 
nunaerical  scheme[l^l  (of  the  second  onler  accuracy  on  the  grid  A)  and  the  numerical  algorithms  are 
described  in  detail  in  the  publications  listed  above.  The  Neumann  boundary  problem  and  the 
ovenelaxation  method  for  computation  of  the  sea  surface  topo^phy  has  been  used.  Ifere  well 
cmly  describe  the  main  features  of  the  basic  nsodel  and  shall  list  the  parameters  used,  without 
details  or  discussion  of  the  specific  calculations.  For  horizontal  velocity  components,  the 
boundary  condition  of  zero  velocity  on  the  botmm  and  on  the  vertical  walls  has  been  u^  for  T&S 
-  dte  isolation  on  all  boundruies;  for  w  -  the  "rigid  lid"  condition  at  the  undisturbed  upper 

boundary,  and  ws  0  on  the  bottom  have  b^n  accepted  as  well: 
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The  wind  stress  has  been  chosen  according  to  the  data  of  Hellennani^^.  The  time  step  has 
been  selected  to  be  ^  min.  and  the  values  of  the  coefficients  of  the  vertical  turbulent  exchange  and 
vertical  heat  and  salt  diffusion  have  been  determined  as  1  cm^/sec.  The  horizontal  turbulent 
diffiision  coefficient  has  been  chosen  to  be  0.0264  h4/3  cm^sec  and  the  horizontal  coefficient  of 
eddy  viscosity  has  been  calculated  using  the  Smagorinsky  fonnulantil. 

The  adaptation  calculation  j^ess  has  been  completed  after  2  or  3  days  of  model  time  (see 
Hgure  1),  i.e.  the  adjustment  time  is  relatively  fast  (wiA  respect  to  the  pre^ction  period).  In  this 
way  the  main  advantages  of  the  adaptive  models  have  been  realized:  using  a  short  integration 
period  good  diagnostic  results  have  been  obtained  (adequate  to  the  experimental  data)  and  at  the 
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same  time  some  of  the  enors  in  the  T-S  fields,  leading  to  meso-  and  small-scale  perturbations, 
have  been  removed.  The  hydrological  adjustment  of  the  fields  results  mainly  in  hydrodynamic 
filtration  of  the  noise  (with  respect  to  the  model)  -  Figures  2,  4.  The  basic  mechanisms  of  the 
adaptation  are  the  "scattering"  and  the  dissipation  of  the  inertial  gravity  waves.  These  aren't 
I^ysical,  but  model  waves.  Tbeir  amplitude  exceeds  the  amplitude  of  the  real  waves,  and  their 
wavelengA  depends  on  the  grid  interval.  These  waves  are  initiated  by  the  "noise"  in  the  input  data 
(by  their  incompatibility  with  the  model),  and  they  are  typical  of  the  models,  based  on  "primitive" 
equations  of  motion,  describing  them  (in  contrast  to  Ac  filtered,  quasigeostrophic  modelll^l). 
Hgure  1  shows  Aat  Ae  oscillations  of  energy  and  enstrophy  disappear  after  2  or  3  days  of 
integration.  If  Ae  mtegration  process  is  continued  funher,  the  fields  are  slowly  smooAed  down 
and  modified  towards  the  model  solution.  There  are  no  significant  difierences  between  Ae  results, 
obtained  after  3  and  after  10  days  of  integration-Ae  Afferences  are  much  smaller  Aan  those 
observed  at  this  stage  of  semiAagnostic  calculations. 

The  process  of  adaptation  removes  Ae  small-scale  structures  (edAes  and  meanders),  and 
causes  undesirable  broadenmg  and  weakening  of  Ae  GBSC.  This  effect  is  much  weaker  Aan  Ae 
influence  of  Ae  traAtional  filtration  procedures,  but  is  neverAeless  observable.  For  some  regions 
of  Ae  GBSC  (for  mstance,  near  Cape  Kaliakra),  Ae  grid  st^  is  too  large,  leading  to  "loss"  of  Ae 
current  jet  Well  list  some  characteristic  features  of  Ae  solution  obtained: 

(1)  Ae  Astributions  obtained  after  Ae  adaptation  are  of  regular  and  "legible"  character,  at 
Ae  same  time  all  large-scale  features  of  Ae  initial  dm  have  been  preserved  (Bgure  2); 

(2)  Ae  adaptation  meAod  Aves  better  results  Aan  Ae  traAtional  filtration  procedures, 
used  for  instance  by  Ae  auAors  oft*^;  Ae  charts  Acre  show  Aat  Ae  GBSC  has  been  significantly 
broadened  and  weakened  by  Ae  filtration  (Figures  3  and  4),  which  is  too  high  a  price  for  Ae 
regularity  obtained 


*  *  * 

Annual  current  fields  at  a  depA  of  1(X)  m  are  shown  m  Figure  2.  They  have  been  obtained 
by  diagnostic  and  adaptive  calculations.  It  can  be  clearly  seen  that  all  basic  structures  of  circulation 
are  more  Astinct  after  Ae  adjustment  process.  The  GBSC,  Ae  basic  dynamic  feature  of  Ae  Black 
Sea,  encircles  Ae  whole  sea  within  one  cyclomc  circular  stream.  It  can  be  considered  as  a  frontal 
current,  which  generates  eddies.  Its  structure,  however,  is  quite  Afferent  from  Ae  common 
concept  of  a  smooA,  symmetrical  stream  wiA  two  gyres  at  Ae  western  and  eastern  parts  of  Ae 
sea.  'The  presence  of  Aese  two  gyres  in  Ae  results  of  model  calculations  is  usually  considered  to 
be  a  sufficient  proof  Aat  Ae  model  is  satisfactorily  describmg  Ae  real  circulation.  ‘The  structure  of 
tlw  currents  is  much  more  complex  -  Ae  stream  meanders,  Ae  gyres  are  elongated  and  change  Aeir 
positions  seasonally.  The  number  of  cores  at  Ae  center  of  Ae  cyclomc  gyres  can  exceed  two.  On 
the  background  of  the  general  cyclonic  rotation  a  series  of  anticyclomc  eddies  develops  close  to  Ae 
shore.  Some  of  Aem  are  present  on  Ae  chans  for  all  four  seasons.  AnoAer  important  feature, 
different  fiom  Ae  common  circulation  representations,  is  Ae  fact  Aat  Ae  current  isn't  always  close 
to  Ae  continental  shelf  slope.  Sometimes  (especially  in  regions  wiA  edAes  and  meanders)  Ae 
current  detaches  from  the  shore.  It's  evident  Aat  the  GBSC  is  subject  to  significant  seasonal 
variations  (see  Appendix),  Aough  its  general  circulation  structure  is  qualitatively  constant.  The 
structure  of  Ae  si^ace  currents  is  preserved  quite  well  down  to  large  depAs,  wiA  decreasmg 
velocity  modulus  only.  Gradually  Ae  GBSC  loses  its  jet  character.  A  significant  change  is 
observ^  at  depAs  below  800  m.  For  instance,  at  Ae  center  of  Ae  sea  at  a  depA  of  ~1000  m,  Ae 
direction  of  Ae  main  stream  is  opposite  to  Aat  at  Ae  upper  layer.  This  strong  (relatively!  -  wiA 
respect  to  tiie  oAer  currents  at  Ae  periphery)  stream  to  Ae  east  is  observed  m  all  seasons  and  for  all 
oAer  levels  to  Ae  bottom,  which  suggests  Aat  it  is  a  stable  feature  of  Ae  circulation  at  large 
depAs.  Hus  result  confirms  Ae  hypothesis  of  Neumannll^i  Aat  at  1(X)0  m  depAs  a  current  exists, 
flowing  in  Ae  direction  opposite  Aat  of  Ae  surface  current 
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Figure  2.  Diagnostic  (a)  and  adaptive  (b)  climatic  currents  at  a  depth  of  100  m. 
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Block  seo  climote  aolinily  2  The  Block  sea  climate  soiinity 


Figure  4.  Climatic  salinity  dletrlbutiona  (?/••)  ati the  sea  surface  for  different 


The  results  described  in  the  present  work  show  that  the  method  proposed  can  be  used  to 
perform  hydrodynamic  analysis  of  the  main  hydrological  fields  in  the  Black  Sea.  Our  experience 
proves  that  the  3D  current  fields  obtained  by  hydrodynamical  adjustment  are  a  reliable  basis  for 
numerical  investigations  of  a  wide  variety  of  applied  problems,  connected  with  sea  pollution 
modelling,  using,  for  instance,  the  direct  and  conju  sated  diffusion  problems  and  obtaining  a 
number  of  integral  characteristics  (balance  approach)!^!^].  The  climate  ^ta  produced  from  a  given 
experiment  should  be  taken  as  a  starting  point  for  obtaining  information  from  the  existing  arrays  of 
hydrological  data.  To  further  develop  and  complement  this  information,  some  requirements  have 
to  be  met  by  the  future  surveys.  It's  important  to  increase  significantly  the  statistical  reliability  of 
the  data  for  large  depths.  The  result  will  be  a  more  reliable  description  of  the  deep  water  dynamics. 
It's  also  necessary  to  provide  information  about  the  wind  stress  and  the  atmospheric  pressure 
above  tire  water  surface  of  the  Black  Sea.  The  data  of  Hellermanf^^l,  which  are  widely  used  in 
numerical  calculations,  are  too  generalized,  and  are  quite  different  from  the  real  wind  characteristics 
of  the  Black  sea.  The  distribution  of  surface  temperature  in  winter  turned  out  to  be  the  "noisiest" 
of  all  fields  used.  This  is  easily  explained  by  the  dynamic  character  of  the  processes,  taking  place 
on  the  siffface  in  this  season,  and  by  the  insitificient  survey  activity  in  winter.  It’s  possible  that  for 
future  investigations,  information  be  obtained  from  satellites  and  frrom  ship  weather  stations  as 
well. 


*  ** 

When  the  results  of  numerical  experiments  are  discussed,  comparisons  of  the  tj^e 
"practically  coincitfing",  "qualitatively  similar",  etc.  are  often  usedt^*  The  quantitative 
evaluations  like  "decrease  by  10  to  20  percent"  are  in^licitly  understood  as  similarity  of  the  results 
as  well.  The  imperfections  of  the  models  used  and  the  lack  of  sufficient  quantitative  (and  often  - 
qualitative  as  w^)  understanding  about  the  investigated  processes  are  the  reason  why  10%  to  20% 
variations  in  Ae  results  cannot  be  regarded  as  sufficiently  reliable.  With  the  future  d^elopment  of 
new  and  more  adequate  models  and  broadening  of  our  knowledge  about  the  investigated 
processes,  the  importance  of  these  conclusions  will  increase  (this  is  the  case  for  example  with 
weather  forecasts). 
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III.  APPENDIX 


CLIMATIC  ATLAS  OF  THE  BLACK  SEA 
SEA  LEVEL  TOPOGRAPHY,  CURRENTS,  SALINITY,  TEMPERATURE 
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It.  Ataliaai  (Madi:  200  arorda) 

The  Black  Sea  is  a  neatly  enclosed  ocean  basin,  exhibiting  many  features  common  with  larger  ocean  basins.  Lacking  an  t^ten 
boundary  and  having  a  limited  exchange  with  sources  of  fresh  ^  salt  water,  this  basin  is  an  ideal  laboratory  for  developing  and 
evaluating  numerical  circulation  models.  The  present  report  describes  one  numerical  model  of  the  Black  Sea,  developed  by  Bulgarian 
and  Russian  scientists.  The  new  r^qjroach  has  the  advantages  of  both  diagnostic  models  (incorpmation  of  experimental  data)  and 
prognostic  models  (producing  hydnxlynamical  adjustment  and  filtered  fields).  Successive  triplication  of  diagnostic  and  prognostic 
models  is  used.  The  temperature  and  salinity  fields  obtained  firmn  observations,  and  currents  obtained  from  diagnostic  niodels,  are 
used  as  the  initial  qiproximation  to  the  prognostic  model.  Judicious  selection  of  an  integration  time  prevents  over-smoothing  of  the 
results  while  preserving  the  stability  of  the  solution. 

Using  this  model,  calculations  have  been  made  at  25  levels  over  a  grid  interval  of  025"  (latitude)  by  OJ".  Input  data  consist 
of  neariy  50,000  observations  taken  over  nearly  100  years,  averaged  over  0.5**  by  0.5*  cells.  Seasonal  fields  of  temperature,  salinity, 
and  velocity  form  the  ouqiut  of  these  experiments.  The  results  provide  the  basis  for  various  hypotheses  that  must  be  tested  using 
future  field  observations  and  more  sr^ihisticated  models. 


17.  PeBinia  AiMtyate  a.  D— crtprerv 

1.  numerical  modeling 

2.  general  drculation 
S.BladtSea 

4.  daui  assimilation 

te  MMaHIlfWOpMl-ElMtMl  TMIM 


•.OOSATIPtohVQroiip 


IS.  Awmasay  atMzimnt 

Approved  for  pubUcalion;  distribution  unlimited. 


ia.  SccurSy  CliM  (Thte  Raport) 
UNCLASSIFIED 

21.  Na  c(  Pagaa 

132 

20.  SoewSy  Ctaao  (Tlite  Pago) 

22Pitoa 

AM  AMSruatfaw  «n  RMw** 


OPTIONAl.  PORH  272(4-77) 
(PomMiV  NnS-35) 
OtpaiiiiMnl  o(  CenmMc* 


